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ABSTRACT 

We propose that magnetic reconnection at the heliopause only occurs where the 
interstellar magnetic field points nearly anti-parallel to the heliospheric field. By us- 
ing large-scale magnetohydrodynamic (MHD) simulations of the heliosphere to provide 
the initial conditions for kinetic simulations of heliopause (HP) reconnection we show 
that the energetic pickup ions downstream from the solar wind termination shock in- 
duce large diamagnetic drifts in the reconnecting plasma and stabilize non-anti-parallel 
reconnection. With this constraint the MHD simulations can show where HP reconnec- 
tion most likely occurs. We also suggest that reconnection triggers the 2-3 kHz radio 
bursts that emanate from near the HP. Requiring the burst locations to coincide with 
the loci of anti-parallel reconnection allows us to determine, for the first time, the vector 
direction of the local interstellar magnetic field. We find it to be oriented towards the 
southern solar magnetic pole. 



1. Introduction 

After crossing the solar wind termination shock in 2004 and 2007, respectively, the Voyager 
1 and 2 spacecraft discovered that the downstream flows remain supersonic with respect to the 

and occurs because most 

(~ 80%) of the solar wind energy upstream of the shock transfers to a non-thermal population 



thermal ions. This configuration had been anticipated (Zank et al.\ 1996a 
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whose presence, although not directly detected, can be inferred from extrapolation of the available 
data ( Decker et aL||2008 Richardson et aZ.||2008 ) . The shocked solar wind can then simultaneously 
exhibit subsonic flow with respect to the energetically dominant component of the plasma and 
supersonic flow with respect to the thermal component. 

Pickup ions likely form the bulk of this non-thermal population. Interstellar neutral atoms, 
unaffected by electromagnetic fields and possessing a velocity characteristic of the solar system's 
motion with respect to the local interstellar medium (LISM), ~ 25km/s, can drift inside the 
termination shock where they encounter the ~ 400km/s solar wind. When these atoms ionize 
they suddenly come under the influence of the solar wind's magnetic field, are "picked up", and 
join the wind's outward flow. The high relative velocity of the neutral atoms becomes an effective 
non-thermal temperature with ksT = mv 2 /2 ~ IkeV. Although not directly observed at the 
termination shock, pickup ions have been detected throughout the inner heliosphere. 

Models suggest that the Voyager spacecraft will remain within the heliosphere for ten to 



twenty years (see Table 1 of Opher et al. (2006)) before finally encountering the HP, the boundary 
between the solar system and local interstellar space. There the interstellar magnetic field abuts 
the heliospheric field in much the same way that, within the solar system, the fields of magnetized 
bodies meet the interplanetary magnetic field at magnetopauses. Magnetic reconnection often 



occurs at such interfaces and has been observed at, among other locations, Earth, Mars (Eastwood 



et al. 2008), and Saturn (Huddleston et al. 1998). By analogy, it is also expected to occur at the 



HP (Fahr et al. 1986). 



However, the high energy content of the pickup ions implies a large value for the heliosheath 
plasma /3 (where /5 = SirnksT / B 2 is the ratio of the thermal and magnetic pressures). Observations 



of reconnection at the terrestrial magnetopause (Scurry et al. 1994) and in the solar wind (Phan 



et al. 2009) suggest that, during high j3 conditions, reconnection only occurs between anti-parallel 
magnetic fields, i.e., when the shear angle between the reconnecting fields is ~ 180°. Particle-in- 
cell (PIC) simulations of reconnection have demonstrated why this occurs: diamagnetic drifts in 



plasmas with a high (3 but shear angles < 180° suppress reconnection (Swisdak et al. 2003). 



In this paper we present PIC simulations showing that diamagnetic stabilization also occurs at 
the HP and limits the possible sites of reconnection to those where the interstellar magnetic field 
lies nearly anti-parallel to the heliospheric field. Large-scale MHD simulations of the heliosphere 
can then determine where HP reconnection should occur. 

We find the locus of anti-parallel reconnection sites to be relatively small, so in situ detection of 
reconnection signatures at the HP by either Voyager spacecraft seems unlikely. However, on several 
occasions since their launch, the spacecraft have detected bursts of radio emission at 2-3 kHz 



believed to be associated with the interaction of interplanetary shocks and the HP ( Gurnett et al 



2003 ) . We argue that reconnection triggered by this interaction drives the electron beams producing 



the radio emission. By combining the locations determined for one set of bursts ( Kurth and Gurnett 



2003) with the requirement that HP reconnection be nearly anti-parallel we can constrain the 
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magnitude and orientation of the local interstellar magnetic field (LISM). In particular, the vector 
direction of the LISM (the difference between B and — B) can be easily determined. This is, to our 
knowledge, the first determination of this parameter. 



2. Heliospheric Model 



As part of our investigation we perform large-scale 3D MHD simulations of the heliosphere. 
The code is based on BATS-R-US, a three-dimensional parallel, adaptive grid code developed by 



the University of Michigan ( Gombosi et al. 1994) and adapted for the outer heliosphere ( Opher 



et al. 2004). 



The Sun is the center of our coordinate system with the x axis parallel to the line through, but 
pointing inward from, the nose of the HP. The z axis lies parallel to the solar rotation axis and the y 
axis completes the right-handed triplet. In order to better describe the different sub-environments 
of the heliosphere such simulations model multiple fluids ( Zank et aZ.||1996b Pogorelov et al.|[2 006); 



we include five hydrogen populations, one ionized and four neutral {Opher et al. 2009a). All four 



neutral populations are described by separate systems of Euler equations with corresponding source 
terms describing neutral-ion charge exchange. The parameters for the inner boundary (located at 



30 AU) were chosen to match those used by Izmodenov et al. (2008): proton density at the inner 
boundary of n = 8.74 x 10 -3 cm -3 , temperature T = 1.087 x 10 5 K, and a Parker spiral magnetic 
field with strength B = 2/iG at the equator. Except where otherwise noted the heliospheric field 
orientation corresponds to that during solar cycle 22 (roughly 1986 to 1996), in which in the x — z 
plane through the nose of the heliopause the field in the northern hemisphere points in the — y 
direction. We denote this as Bsw,y < 0. 

The outer boundary conditions are n = 0.06 cm -3 , velocity = 26.3 km/s, and T = 6519 K. 
The neutral hydrogen in the LISM is assumed to have n = 0.18 cm -3 and the same velocity and 
temperature as the ionized LISM. We use fixed inner boundary conditions for the ionized fluid and 
soft boundaries for the neutral fluids. We impose outflow outer boundary conditions everywhere 
except the — x boundary where inflow conditions are used for the ionized and neutral populations 
from the interstellar medium. The grid's outer boundaries are at ±1500 AU in all three directions 
and the computational cell sizes range from 0.73 to 93.7 AU. The orientation of the interstellar 
magnetic field is characterized by the angles ais, the angle between the field and the interstellar 
wind, and j3is, the the angle between the field and the solar equator. Given ajs and fiis, the vector 
direction of the field is still undetermined up to a sign that can be fixed by noting whether the y 
component of &ism is positive or negative. 

In Figure [T] we show an overview of a heliospheric simulation with BjsM,y > and the he- 
liospheric polarity of solar cycle 22. The colors represent |B| and the black lines are streamlines 
of the flow. The roughly circular surface at a radius of 80 AU from the Sun where the magnetic 
field strength increases from ~ 0.02nT to ~ 0.2nT is the termination shock. Downstream from 
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the shock the plasma density and temperature increase and the solar wind velocity decreases, as 
expected. The neutral hydrogen species do not have a dramatic effect on the overall morphology 
of the system, although they affect the distances to the termination shock and HP, bringing them 
closer to the Sun, as can be seen by comparison to previous simulations employing only the ionized 



fluid (Opher et al. 2006, 2007). 



Reconnection occurs at the grid scale at a rate determined by numerical details rather than 
by physical processes in these MHD simulations. Previous modeling has shown that without the 
addition of some sort of extra-physical diffusive process, e.g., an anomalous resistivity, MHD recon- 
nection is inherently slow ( Biskampf l986 ) . The inclusion of the Hall term in the generalized Ohm's 



law dramatically enhances the reconnection rate {Birn et al. 2001), but necessitates the resolving 



of length scales on the order of the ion inertial length di = c/uj p i. For a typical heliopause density 
of 0.05 cm" 3 , di ~ 10 8 cm ~ 10~ 5 AU, several orders of magnitude smaller than the resolution of 
the global simulation. Since the MHD model does not correctly describe reconnection, we use it as 
a starting point for kinetic simulations of HP reconnection. 



Kinetic Reconnection Simulations 
3.1. Numerical Methods 



For our kinetic simulations we use p3d, a massively parallel PIC code ( Zeiler et aL||200"2" ). We 



use a different coordinate system than the MHD simulations because our computational domains 
represent only a tiny, arbitrarily oriented, portion of the MHD domain. The inflow into and outflow 
from an X-line lie parallel to y and x, respectively, while the reconnection electron field and any 
magnetic guide field parallel z. In the simulations presented here we assume out-of-plane derivatives 
vanish, i.e., d/dz = 0. Although this choice eliminates any structure in the z direction, previous 



studies indicate the basic features of reconnection remain unchanged (Hesse et al. 2001). 



Masses are normalized to the ion mass m,, magnetic fields to the asymptotic value of the 
reversed field -Bo, and the density to the value at the center of the current sheet. Other normal- 
izations derive from these: velocities to the Alfven speed va, lengths to the ion inertial length 
di = c/ojpi (with ujpi the ion plasma frequency), times to the inverse ion cyclotron frequency fi" 1 , 
and temperatures to m,iV 2 A . 

We set the speed of light (in normalized units) to 20 and the ratio of the ion and electron masses 
to m,i/m e = 100. The spatial resolution is such that there are > 4 gridpoints per electron inertial 
length and ~2 per Debye length. The Courant condition determines the particle timestep and we 
substep the advancement of the electromagnetic fields. A typical cell contains ~ 100 particles and 
our simulations follow > 10 9 particles. All of our runs conserve energy to better than 1 part in 100. 
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3.2. Simulation Results 

The components of any plasma that includes both a magnetic field and a non-parallel pressure 
gradient undergo a diamagnetic drift given by 



Vpj X B 



(1) 



where pj = njksTj is the pressure and qj is the charge of species j. Note that, because of their 
charges, ions and electrons drift in opposite directions and that separate populations of the same 
species (in particular, pickup and thermal ions) can drift at different speeds. In our case the 
pressure gradient normal to the HP crossed with the guide magnetic field produces drifts parallel 
to the reconnection outflows. 



Using PIC simulations Swisdak et al. (2003) demonstrated that under such conditions the 



X-line convects in the ion rest frame with a speed given by |v* e | + |v*j|, the sum of the electron 
and ion diamagnetic drifts. If the drift velocity of the X-line exceeds the speed of the nominally 
Alfvenic outflows, reconnection is suppressed. Qualitatively, suppression occurs when the X-line 
propagates fast enough that the plasma does not have sufficient time to establish the necessary 



flow configuration for reconnection before the X-line passes. Quantitatively, Swisdak et al. (2003) 
proposed that diamagnetic drifts suppress reconnection when 

v*,j > vaj (2) 

where vaj is the Alfven speed for a particular populations (i.e., calculated based on nj) corre- 
sponding to the reconnecting magnetic field. This condition can be derived from balancing the 
centrifugal and magnetic tension forces exerted on a fluid streaming with velocity v#j towards the 
X-line (see the Appendix for details). When is sufficiently large the magnetic tension cannot 
eject the field line from the vicinity of the X-line and reconnection stalls. 

Equation [2] can be reformulated as a condition relating the jump in the plasma parameter 
across the current layer, A/3, and the shear angle 9 between the reconnecting fields: 

21 

Afo > --^tan(0/2) (3) 

where L p represents a typical pressure scale length near the X-line. Both simulations and observa- 
tions at the Earth's magnetopause suggest that L p /di ~ 0(1) ( Berchem and Russell\1982 Eastman 



et al. 1996). According to equation pi only anti-parallel reconnection, 9 tt, occurs for /5 > 1. 



Observations of reconnection at the magnetopause { Scurry et aZ.||1994 ) and in the solar wind (Phan 



et al. 2009) support this conclusion. 



The inferred large energy content of pickup ions in the heliosheath implies a substantial value 
for /3, which, when coupled with equation [3j leads to the conclusion that only anti-parallel recon- 
nection occurs at the HP. To check, we performed PIC simulations of HP reconnection at two shear 
angles. Cuts through the HP at two places in the simulation of Figure [T] established the initial 
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equilibria for these simulations. Due to the relatively large grid size the heliopause in the global 
simulations has a thickness of approximately 10 AU. This is only an upper limit, however, as the 
real thickness is certainly closer to the ldi <C 10AU observed at the terrestrial magnetopause. Fur- 
thermore, in both experiments ( Yamada\2007 ) and simulations ( Cassak et al.|2005 ) reconnection is 
seen to proceed slowly in thick current sheets before accelerating rapidly when the sheet thickness 
reaches d%. In our kinetic simulations the MHD models provide the asymptotic parameters of the 
reconnecting plasmas, but we take the width of the current layer to be ldi. In Figure [2] we show the 
initial profiles of the densities and temperatures at the two locations used to provide the asymptotic 
profiles for the PIC simulations. 

The initial proton distribution is a superposition of a cold, Maxwellian population representing 
the solar wind and a much hotter Maxwellian (20% by number) corresponding to the pickup ions. 
Choosing a thermal distribution for the pickup ion distribution is a somewhat crude approximation. 
The actual distribution is still a subject of investigation since neither of the Voyager spacecraft can 
measure it directly. Close to the Sun, where many of the pickup ions are created, a ring distribution 
could be appropriate, but as the solar wind expands outward a number of altering effects have been 



identified (Isenberg 1987, Chalov et al. 1995) even before the interaction with the termination 



shock introduces further modifications (Zank et al. 2010). In our kinetic simulations the particle 



distribution functions are free to evolve (and do so, particularly in the reconnecting current layer) 
and so the effects of choosing a different initial distribution are ultimately unclear, albeit of obvious 
interest. 

The first site has nearly anti-parallel reconnecting fields (shear angle 9 = 165°), while in 
location 2 the fields have a shear angle of 9 ~ 130°. In order to specify the plane of reconnection 



for these 2D simulations we use the prescription derived in Swisdak and Drake (2007) for finding 



the X-line orientation based on the asymptotic plasma parameters. At locations 1 and 2 the 
values of A/3 are ~ 2 and 4, respectively. Equation [3] then predicts that reconnection should be 
diamagnetically stabilized at location 2, but not at location 1. 

In Figure [3^, we show a snapshot of the out-of-plane current density from the anti-parallel 
case. The bulge of the separatrices into the low field strength heliosheath is typical of asymmetric 



reconnection and has been observed in simulations of the terrestrial magnetopause ( Krauss- Varban 



et a/. 1 1999 ; Nakamura and Scholerf2000 ) . Since the solar wind has the lower field strength in both 
the terrestrial and heliospheric systems the separatrices bulge toward the Sun in both cases. 

We show a similar image from the non-anti-parallel reconnection case in Figure [3j3. Because of 
the diamagnetic drift the X-line has drifted to the right from its initial position at approximately 
the velocity given by equation [TJ A comparison of the size of the islands downstream from the 
X-lines in the two simulations indicates that essentially no flux has reconnected in Figure [3jx The 
normalized reconnection rate as a function of time for the two simulations can be seen in Figure 
[3fc. The anti-parallel case (solid line) clearly exhibits reconnection at a rate of ~ 0.06, typical of 



"fast" Hall reconnection (Birn et al. 2001). In the second case (dashed line), where reconnection 
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is not anti-parallel, very little flux reconnects. 



4. Radio Emission 



The plasma wave instruments on the Voyager spacecraft have detected several bursts of radio 
emission at frequencies between 2 and 3 kHz during their more than 25 years of operation. A strong 
interplanetary shock passed by the spacecraft before each burst, which led to the suggestion that 
the interaction of global merged interaction regions (GMIRs) and the HP boundary produced the 



emission ( Gurnett et al. 2006 ) . An application of several direction-finding techniques to one set of 



bursts demonstrated that the sources lie at the radial distance of the HP and stretch along a line 
parallel to the galactic plane and passing near the nose of the heliosphere. 



Cairns and Zank (2002) have suggested a detailed mechanism for producing the radio emission. 



Just beyond the HP, lower hybrid waves can resonantly accelerate ambient electrons that then 
serve as a seed population for shock acceleration. The enhanced electron beams that form then 
excite Langmuir oscillations, which subsequently decay, via nonlinear effects, into the observed 
electromagnetic radiation. This scenario successfully explains, first, why the Voyager spacecraft do 
not observe radio emissions (or the precursor Langmuir waves and electron beams) when passed by 
shocks in the inner heliosphere, and second, why the 2-3 kHz signal seems to be radially confined 
to a thin layer near the HP. 

Reconnection at the HP can also explain these observations. Within the heliosphere the current 
layers composing the heliospheric current sheet are broad enough (~ lOOdj) that, despite occasional 
observations of X-lines (Gosling 2007), reconnection proceeds slowly on average. Instead, the 



most significant reconnection should occur when the plasma encounters a boundary, as happens at 
planetary magnetopauses and, presumably, the HP. The superposition of solar wind structures in 
a GMIR means that the magnetic field strength can be larger, by a factor of ~ 4, than the ambient 
value ( Whang et al.\19 95). At the heliopause this enhanced field will produce stronger reconnection 
and therefore naturally accounts for the association between the arrival of the GMIRs at the 
heliopause and the bursts of emission. Second, observations and simulations have demonstrated 
that reconnection can produce both beams of electrons traveling at the electron Alfven speed 



( Cattell et al. 2005 ) and more isotropic populations at relativistic energies ( 0ieroset et al. 2002 



Hoshino et al.|200T Drake et a/.|2006 ) . The electromagnetic decay mechanisms described in Cairns 



and Zank (2002) can then generate the 2-3 kHz radiation from these electrons. Reconnection, then, 



naturally localizes the emission at the heliopause, produces the necessary energetic electrons, and 
correlates the emission with the arrival of GMIRs. 

A significant test of this model is whether the locations of the 2-3 kHz radiation, as determined 



by \Kurth and Gurnett (2003) for a series of events between 1992 and 1994 (during solar cycle 
22), correspond to the regions where reconnection should occur. The diamagnetic stabilization 
of reconnection facilitated by pickup ions described in Section [3] implies that reconnection will 
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only occur where the reconnecting fields are anti-parallel. The locations where such reconnection 
occurs depend on the orientation of the interstellar magnetic field, which will not be directly 
measured until one of the Voyager spacecraft passes the heliopause. However inferences from 
indirect measurements, including backscattered Lyman-a emission and the difference in the flow 
directions of interstellar hydrogen and helium (the hydrogen deflection plane) do present some 



constraints (see Opher et al. (2007) and references therein). Other constraints have been derived 



based on the heliospheric asymmetries and heliosheath flows measured by Voyager 2 ( Opher et al. 



2009b). If reconnection is the source of the 2-3 kHz radiation, the requirement that the sources lie 
near the locus of anti-parallel reconnection can provide a strong further constraint on the vector 
direction of the interstellar magnetic field. 

We used the results of the global MHD model described in Section [2] to compare the locations 



of the radio sources determined by Kurth and Gurnett (2003) and the locations of anti-parallel 



reconnection at the HP. We first define the HP in the simulation as the empirically defined surface 
where log(T) = 10.9; other values shift the HP somewhat but do not significantly affect our 
conclusions. Assuming the field in the solar wind, T$sw is purely azimuthal upstream of the HP, 
and hence neglecting any draping of the heliospheric field, we calculate where T$sw ' (&ISM X 
n)/|B/sjy/||B,svH vanishes, where Bism is the field in the local interstellar medium and n is the 
normal to the HP. (We use this metric rather than the simpler (Bism''Bsw)/\'Bism\\'Bsw\ = ~~ 1 to 
eliminate the effects of spurious normal magnetic fields at the heliopause.) This quantity vanishes 
for both anti-parallel and parallel field configurations so we must combine it with another measure 
(we use the field shear angle) to identify the anti-parallel locations. 

In Figure [4^, we show the map produced by the interstellar field parameters — ais = 30°, 
Pis = 60°, |B| = 4.4 fiG, Bism,v > — that provide the best fit. These parameters are consistent 
with the values determined by other methods. Our models suggest that the locus of anti-parallel 
reconnection is more sensitive to a is and |B| than it is to Pis- Changing ais by 10° or |B| by 
20% moves the anti-parallel locus significantly away from the sources, while shifts in Pis of up to 
20° keep the concordance acceptable. Panels and|4j; displays the map for ais = 30°, Pis = 80° 
and ais = 20°, Pis = 60°, parameters at the edge of what we deem the acceptable range. 

Perhaps of most interest is the dependence of our results on the vector sign of the local 
interstellar magnetic field. Figure [4] corresponds to an interstellar field with BisM,y > 0. Reversing 
the sign of ~Bism leads to a heliopause with no anti-parallel reconnection sites near the nose. If 
anti-parallel reconnection^] causes the 2-3 kHz radio bursts at the heliopause, then the existence of 
these bursts fixes the vector direction of the interstellar magnetic field. 

The polarity of the heliospheric field reverses due to the 11-year solar cycle. A location that 
exhibits anti-parallel reconnection during one cycle will, in the next 11 years, find that the he- 



1 Although diamagnetic stabilization provides one way, any mechanism that prevents non-anti-parallel reconnection 
will have the same effect. 
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liospheric and interstellar fields are parallel, and hence unlikely to reconnect)^} A similar change 
happens at the Earth's magnetosphere when the interplanetary magnetic field rotates (on a much 
faster time scale) from, for instance, southward to northward. However, just as in the magneto- 
spheric case, even when anti-parallel reconnection is not possible near the nose of the HP it can 
occur, albeit more weakly, at much higher latitudes near the cusps of the Earth's field. This effect 
cannot be directly tested with the data from Figure [4] because all of the source locations were 
determined during a two-year period (1992-1994) of solar cycle 22. The Voyager spacecraft did 
detect radio bursts during solar cycles 21 and 23, but could not determine the locations of their 
sources, specifically if they occurred at the cusps as our model suggests. Intriguingly, these events 
were substantially weaker than those of solar cycle 22, which agrees with the observation that in 
the terrestrial case cusp reconnection is weaker than its equatorial counterpart. 



Figure [5] displays results for the solar cycle 23 polarity of the solar dipole (Bsw,y > 0) for 
a case otherwise identical to that shown in Figure U — /3js = 60°, ajs = 20°, |B| = 4.4//G, 
BiSM,y > 0. If plotted in the same manner as Figure [4] the entire figure would be white because the 
interstellar and heliospheric fields have shear angles < 90°. Instead the figure shows that T$sw an d 
B/sm are nearly parallel over the entire face of the heliopause. In this case we expect anti-parallel 
reconnection to only happen at high heliospheric latitudes, however the inadequate resolution of 
the global models there does not allow us to accurately predict the favorable locations. 

Future data sources may allow us to test this dependence. This model predicts a transition 
to energetic reconnection near the nose in solar cycle 24 as the polarity of the heliospheric current 
sheet reverses again. However current observations suggest that the Sun may be entering a deep 
solar minimum, with very little magnetic activity. If true it may be accompanied by a concomitant 
decrease in heliopause reconnection. In any case, in situ samplings of reconnection by the Voyager 
spacecraft, aside from being unlikely due to the small odds that a given outward trajectory through 
the heliopause will pass near an X-line, will not be possible for at least another decade. Indirect 



measurements of reconnection at the HP are more likely in the near future. Observations (Lin 



et aZ.||2003 ) and simulations ( Drake et aZ.||2009 ) have demonstrated that reconnection can energize 



ions (this occurs in addition to the electric energization necessary for the 2-3 kHz radio emission). 
When these energetic ions interact with the surrounding plasma they can undergo charge exchange 
and create energetic neutrals that can perhaps be sensed remotely by such missions as IBEX. Such 
a signal would appear in addition to the band of energetic neutrals recently reported by IBEX 
\McComas et q/.||2009| > which is, in this view, unrelated to HP reconnection. 
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in the paper. Computations were carried out at the NASA AMES Research Center and National 



2 Note that in the MHD code reversing the orientations of both the interstellar and heliospheric fields has no 
effect on the loci of anti-parallel reconnection: B v jsm > with a given orientation of the heliospheric field behaves 
identically to B y jsM < with the opposite heliospheric orientation. 



Energy Research Supercomputing Center. M. O. and F. A. B. acknowledge the support of NASA- 
Voyager Guest Investigator grant NNX07AH20G and the National Science Foundation CAREER 
grant ATM-0747654. 



A. Condition for Diamagnetic Suppression of Reconnection 

To more formally derive the condition for diamagnetic suppression of reconnection given in 
equation [2j begin by considering the outflow from an X-line. In the simplest case, with no diamag- 
netic drifts, reconnection produces bent magnetic field lines that accelerate away from the X-line 



due to the J X B force. Specifically, in the coordinate system described in Section [3TT| where outflow 
is parallel to ±x, the force is proportional to J z B y . 

In the case with diamagnetic drifts along the x axis one of the outflows from the X-line will 
be in the direction opposite to the drifting plasma. The change in momentum that a bent field line 
can cause in a time At over a box of dimensions A x by A y is 

J z B y A x A y At/c (Al) 

The unbending of the field line, and hence reconnection, can only occur if this quantity is large 
enough to overcome the momentum of the plasma in the box traveling with the diamagnetic velocity 

— pjV*jA y (v*jAt) (A2) 

where pj and v*j are the density and diamagnetic velocity of species j. Combining these expression 
leads to a condition for the suppression of reconnection. By using Ampere's law to substitute 
for J z = (dB x /dy)/ (Air/c) and Gauss's law to equate dB x /dx and —dB y /dy we arrive at the 
suppression condition: 



B 2 



Pi<i > t (A3) 
Equation A~3 is equivalent to equation [2] of the main text. 
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X AU 

Fig. 1. — Cut in the x — z plane of a 3D MHD simulation of the outer heliosphere with BisM,y > 
and a heliospheric field with the solar cycle 22 polarity (Bsw,y < 0). Colors denote magnetic field 
strength (in nT) and the black lines represent flow streamlines. The HP traces the outer edge of 
the the heliospheric current sheet, which is shown by the blue line paralleling the — x axis and then 
deflecting northward in the heliosheath. 
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Fig. 2. — Cuts taken through the HP from the MHD simulation showing the density and pickup 
ion temperature. Panels (a) and (c) come from a site, location 1, with anti-parallel reconnection; 
panels (b) and (d) from location 2 where the fields have a smaller shear angle. 
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Fig. 3. — Panels (a) and (b): Out-of-plane current density overlaid by magnetic field lines for two 
PIC reconnection simulations. The heliosheath and LISM plasma are above and below the current 
sheet respectively. Panel (a) corresponds to reconnection at location 1, where the fields are anti- 
parallel. Note the well-developed X-line. In panel (b) we show reconnection at location 2, where 
the fields are not anti-parallel. The X-line, initially at x/di = 0, has drifted due to diamagnetic 
effects. Panel (c): Reconnected flux versus time for the simulations shown in panels (a) (solid line) 
and (b) (dashed line). The reconnection rate is given by the slope of each line. 
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Fig. 4. — Locations of anti-parallel reconnection at the HP for three configurations of the 
interstellar magnetic field as seen from outside the heliosphere looking inward. SELat and 
SELong are solar ecliptic latitude and longitude, respectively. Colors denote the quantity 
E$5ty • (Bjsm X n)/|B/5A/||B5'vi / | and have been saturated at either extreme. White regions are 
places where the shear angle between J5sw and &ism is less than 90°. The squares represent the 
sources of radio emission, as determined by Kurth and Gurnett (2003). In panel (a) j3is = 60° and 
ais = 30°, in panel (b) fiis = 80° and ais = 30°, and in panel (c) fiis = 60° and ais = 20°. In 
all panels BjsM,y > and B$w,y < (solar cycle 22 polarity). 
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Fig. 5. — Angles between the interstellar and heliospheric fields for the solar cycle 23 orientation 
of the solar dipole. As in Figure [4] colors denote T$sw ' (B/SAf x n)/|B/s-jv/| iB^wl , but here the 
fields have shear angles < 90°. For this run, /3js = 60°, ajs = 20°. 



